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ABSTRACT 

We studied and compared the radial profiles of globular clusters and of the stellar 
bulge component in three galaxies of the Fornax cluster observed with the WFPC2 of 
the HST. A careful comparative analysis of these distributions confirms that stars are 
more concentrated toward the galactic centres than globular clusters, in agreement 
with what was observed in many other galaxies. If the observed difference is the result 
of evolution of the globular cluster system sstarting from initial profiles similar to 
those of the halo-bulge stellar components, a relevant fraction of their mass (74%, 
47%, 52% for NGC 1379, NGC 1399 and NGC 1404, respectively) is disappeared in 
the inner regions, likely contributing to the nuclear field population, local dynamics 
and high energy phenomena in the primeval life of the galaxy. An indication in favour 
of the evolutionary interpretation of the difference between the globular cluster system 
and stellar bulge radial profiles is given by the positive correlation we found between 
the value of the mass lost from the GCS and the central galactic black hole mass in 
the set of seven galaxies for which these data are available. 

Key words: galaxies: globular clusters; galaxies: star distribution; galaxies: globular 
cluster evolution; galaxies: nuclei 



1 INTRODUCTION 

A lot of elliptical galaxies have globular cluster systems 
(hereafter GCSs) that are less concentrated to the galactic 
centre than the bulge-halo stars, thus showing radial distri- 
butions with larger core radii than those of the bulges. Two 
classical examples of this phenomenon are the two Virgo 
giants M87 and M49 (Grillmair et al., 1986, McLaughlin, 
1995), to which our Galaxy and M31 have been added later 
(Capuzzo-Dolcetta & Vignola, 1997). The sample of galaxies 
that shows this particular feature has been recently enlarged 
thanks to 11 elliptical galaxies selected by Capuzzo-Dolcetta 
& Tesseri (1999) in a sample of 14 galaxies with kinemat- 
ically distinct cores (Forbes et al.,1996) observed with the 
WFPC2 of the Hubble Space Telescope. Though we could 
not yet conclude that this phenomenon is common to all 
galaxies, we may say that there is no case where the halo 
stars are less concentrated than the GCS. Moreover there is 
a general agreement on that the difference between the two 
radial distributions is real and not due to a selective bias. 
Consequently, different hypotheses have been advanced with 
the purpose to explain this feature. Among these, two seem 
the most probable: 

(i) the difference between the two distributions reflects 
different formation ages of the two systems, as suggested by 



Harris & Racine (1979) and Racine (1991); in their opinion 
globular clusters originated earlier, when the density dis- 
tribution was less peaked. However, this hypothesis cannot 
explain why the two distributions are very similar in the 
outer galactic regions. 

(ii) another explanation is based on the simple and rea- 
sonable assumption of coeval birth of globular clusters and 
halo stars (McLaughlin, 1995), with a further evolution 
of the GCS radial distribution, while the collisionless halo 
stands almost unchanged. 

The GCS evolution is due to 

(i) dynamical friction, that brings massive clusters very 
close to galactic centre, 

(ii) tidal interaction with a compact nucleus (see for ex- 
ample Capuzzo-Dolcetta, 1993). 

The combined effect of these dinamical mechanisms acts to 
deplete the GCSs in the central, denser regions, leaving al- 
most unchanged the outer profile, that so remains similar 
to the halo stars' one. The efficiency of the mentioned phe- 
nomena is higher in galactic triaxial potentials (see for exam- 
ple Ostriker et al., 1989, Pesce, Capuzzo-Dolcetta & Vietri, 
1992), in which we find a family of orbits, the 'box-like', that 
do not conserve any component of their angular momentum 
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and then well sample the central galactic regions. Pesce et al. 
(1992) showed as globular clusters moving on box orbits lose 
their orbital energy at a rate an order of magnitude larger 
than on loop orbits of comparable size and energy. These 
results showed as the previous evaluations of the dynamical 
friction efficiency, based on very simplified hypotheses on 
the globular cluster orbit distribution, led to an overstimate 
of the dynamical braking time scales. On the contrary, it 
has been ascertained that massive globular clusters in tri- 
axial potentials could be braked during their motion so to 
reach the inner regions in a relatively short time(Capuzzo- 
Dolcetta, 1993). There they could have started a process of 
merging, giving origin to a central massive nucleus (eventu- 
ally a black hole) or fed a preexistent one (Ostriker, Binney 
& Saha, 1989; Capuzzo-Dolcetta, 1993). Under the hypoth- 
esis that the initial GCS and halo-bulge radial distributions 
were the same, an accurate analysis of the observations al- 
lows an estimate of the number of "missing clusters" and 
therefore of the mass removed from the GCSs. 

Actually, McLaughlin (1995), Capuzzo-Dolcetta & Vig- 
nola (1997) and Capuzzo-Dolcetta & Tesseri (1999), scaling 
the radial surface profiles of the halo stars of a galaxy to that 
its globular cluster system, estimated the number of missing 
globular clusters as the integral of the difference between 
the two radial profiles. A study like this led, for example, 
Capuzzo-Dolcetta & Vignola (1997) and Capuzzo-Dolcetta 
& Tesseri (1999) to suggest that the compact nuclei in our 
galaxy, M31 and M87 as well as in many other galaxies could 
have reasonably sucked a lot of decayed globular clusters in 
the first few Gyrs of life. The aim of this paper is an ex- 
tension of previous (McLaughlin 1995; Capuzzo-Dolcetta & 
Vignola 1997; Capuzzo-Dolcetta & Tesseri, 1999) works to 
the globular cluster data of a sample of three elliptical galax- 
ies in the Fornax cluster (NGC 1379, NGC 1399 and NGC 
1404) whose HST data were taken by Forbes et al. 1998a,b, 
hereafter referred to F98a and F98b. 



2 THE DATA AND THE RESULTS 

In this Section we present and discuss our analysis of GCS 
radial distribution in the overmentioned galaxies in the For- 
nax cluster. The data come from recent HST WFPC2 obser- 
vations (F98a, F98b) and, so, they take the advantage of the 
HST excellent resolution, which allows an accurate analysis 
of the crowded inner regions of galaxies at Virgo and Fornax 
cluster distance. Moreover it allows the subtraction of most 
of background galaxies. The available density data are, as 
we will see later, at a high level of completeness. 
We have approximated the observed density profiles with a 
modified "core model" 
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where Eo, 7 and r c are free parameters (r c is called core 
radius, but, to be precise, it corresponds to the usual defini- 
ton of being the distance from the centre where the surface 
distribution halves its central value just when 7 = 1. We 
have chosen the free parameters in (1) minimizing the r.m.s 
error, i.e. 



where fi e /o; are, respectively, the fitted and the observed 
values of the surface density. We have defined the weights pi 
as the inverse of the error bars of the surface density data. 
The modified core model approximates very well the GCS 
density data in the inner regions of the two galaxies NGC 
1399 and NGC 1404 (a ~ 10~ 4 ), while the a value results 
a little higher (~ 10~ 3 ) for NGC 1379. As we will see in 
the following, this is justified by the somewhat peculiar be- 
haviour of the light profile of the latter galaxy. By means 
of fitting formula (1) we estimated the number of "missing" 
clusters in a galaxy and the corresponding mass removed 
from its GCS. In fact, an approximated value of the mass 
fallen to the galactic centre is obtained through the two val- 
ues, Ni and {mi), of the number and of the mean mass of 
the missing ('lost ') globular clusters. A priori, the deter- 
mination of (mi) needs a detailed evaluation of the effects of 
evolutionary processes on the original mass spectrum of the 
GCS. However, the most relevant phenomena (tidal shocking 
and dynamical friction) act on opposite sides of the initial 
mass function then - if the initial mass function is not too 
asymmetric - we expect that the mean value of the glob- 
ular cluster mass has not changed very much in time (see 
Capuzzo-Dolcetta & Tesseri, 1997) . Hence we can assume 
the present mean value of the mass of globular clusters, (m) , 
as a good reference value for (mi). To estimate (m) it is thus 
necessary to assume a mass function that satisfactorily rep- 
resents the present distributions of globular clusters in the 
three galaxies of our sample. On the basis of the observed 
similarity among the mass spectra of elliptical galaxies (Har- 
ris & Pudritz 1994; McLaughlin 1994), we have obtained a 
single mass spectrum for the Fornax galaxies normalizing 
the mass spectrum adopted by McLaughlin (1995) for the 
Virgo giant M87. We have made this normalization using a 
mass to light ratio (3 = (M/L) VQ equal to 1.5 for the GCS. 
The resulting mass spectrum is 



N(m) oc m 
with s defined as 

r s = 0.5 ± 0.2 
s=< si = 1.65 ±0.10 
I s 2 = 3.0±0.3 



(3) 



m min < m < 1.2 x 10 5 M Q 
x 10 5 M X < m < 1.5 x 10 6 j 
1.5 x 1O 6 M < m < m max 



where m min and m max depend explicitly on /?.Of course an 
error in the value of s leads to an error in the evaluation of 
the mean mass (m). Then we think important to check the 
degree of reliability of our method by estimating (in the way 
described in the Appendix A) the mean error made using the 
mass spectrum (3) to evaluate the mass. 



2.1 NGC 1379 

NGC 1379 is an EO galaxy in the Fornax cluster studied by, 
e.g., Harris & Hanes (1987), Kissler-Patig et al. (1997) and 
later F98b. Harris & Hanes (1987) compared the radial pro- 
file of the GCS in NGC 1379 with the surface brightness of 
the galaxy itself as given by Schombert (1986). Their study 
was limited to the radial range 5-35 kpc, and they did not 
detect any difference between the two profiles. This was con- 
firmed by Kissler-Patig et al. (1997) that used ground-based 
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data over the 3-10 kpc range. The recent analysis of F98b, 
based on data of the HST WFPC2, confirms the similarity 
between the two surface profiles in the outer regions but, at 
the same time, show as the two differ in the inner regions. 
Thanks to these latter data, that cover a radial range larger 
than previously available, we managed to make a compar- 
ative study between the halo stars' surface profile and the 
GCS one. 

2.1.1 Observed data and the modified core model 

In F98b detailed plots of the radial profile of the GCS of 
NGC 1379 are presented, thanks to the high resolution of 
HST. The profiles shown in the paper are two: one (contain- 
ing GCs with B < 25.5) is complete and uncontaminated and 
the other (composed by objects with B> 26.5 and without 
correction for incompleteness) is expected to be composed 
mostly by background galaxies. 

The radial profile of the brighter sample begins to de- 
crease from ~ 10 to 80 arcsec (~ 1 — 7 kpc); at 80 arcsec 
from the centre the profile is lost in the background (see 
F98b). 

Fig.l shows the radial profile of clusters with B < 25.5 
with a background of 2.7 objects per arcmin 2 subtracted. 
The superposed dotted line represents the luminosity profile 
of the underlying galaxy as provided by Kissler-Patig et al. 
(1997), scaled vertically to match the GCS profile in the 
outer regions. In this way the two profiles agree well at ~35- 
70 arcsec (3-6 kpc). The logarithmic slope of the GCS profile 
is ~ —2.4 at f>35 arcsec. Inwards of 30 arcsec it flattens 
out. 

F98b estimated a core radius of 23±6 arcsec (2.0±0.5 
kpc) assuming it as the distance from the centre where 
the surface density halves its central value. In its turn, the 
galaxy light seems to decline with a constant slope from ~ 10 
to ~ 50 arcsec, where it seems to change slope slightly. To 
be sure of completeness, we have initially analysed only the 
brighter sample (B< 25.5): in a second analysis, we added 
the clusters of the fainter sample. Fitting the GCS density 
data with the law (1), we have the best fit with Eo = 220.62 
arcmin -2 , 7 = 1 and r c = 23 arcsec (solid line in Fig.l). The 
modified core model approximates very well the density data 
in the outer regions, while it seems to be an overestimate in 
the inner regions. According to F98b the central dip of the 
observed GCS radial distribution is not an artefact of their 
analysis, but it represents a real drop in the volume density 
of clusters near the nucleus of galaxy. In fact they estimated 
a 97% completeness of the density data in the galactic region 
inwards of ~10 arcsec. 

The observed decreasing behaviour toward the galactic 
centre of the GCS distribution of NGC 1379 seems to indi- 
cate that the dynamical depleting mechanisms are efficient, 
in this galaxy, up to significantly large radii. 

2.1.2 The "missing clusters" in NGC 1379 

As we said in the Introduction, the estimate of the (poten- 
tially) missing clusters can be obtained as the integral of 
the difference between the GCS and halo stars' surface den- 
sity distribution, after this latter has been normalized. For 
this reason we have considered only the radial range where 
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Figure 1. Surface density profile for the the sample of globular 
cluster system in NGC 1379 (black dots) with B < 25.5 from 
F98b. The solid line is the modified core model fit for GCS made 
in our study. The dotted curve is the surface brightness profile of 
the underlying galaxy, arbitrarly normalized to match the radial 
profile of cluster system in the outer regions. Radial distance is 
in arcsec (1 arcsec = 89 pc); surface density in arcmin" 2 . 

the two profiles differ, i.e. 1.77-37.58 arcsec (see Fig.l). The 
bulge-halo radial profile has an almost constant slope, -2.4, 
in the range 10 to 50 arcsec from the centre, but for our 
purpose it has been necessary to know the star profile also 
in the region within 10 arcsec. This has been made possible 
thanks to Grillmair, who gave us (1999, private communi- 
cation) the model used to subtract the stellar light from 
globular clusters profile (99 % complete for B < 25.5). We 
obtained a logarithmic slope of the stellar profile of about -1 
in the radial range 1.77-10 arcsec (see Fig.l). In the annulus 
1.77-37.58 arcsec, the number of stars is 512, on the other 
side, the modified core model fitting the globular cluster 
distribution gives 132 as number of clusters observed in the 
same annulus. Then, in the assumption that the normalized 
stellar profile gives the initial distribution of globular clus- 
ters, we find that the number of "missing clusters" is 380. 
Hence NGC 1379 could have lost about 74% of its initial 
population of globular clusters. 

2.1.3 Data completeness and the mass fallen to the 
galactic centre 

The relevant number of clusters possibly missing in the inner 
region of NGC 1379 corresponds to a great amount of mass 
supposedly fallen close to the galactic centre. Remembering 
that the value obtained of Ni is restricted to clusters with 
B < 25.5, we decided to provide two different evaluations of 
the mass removed from the GCS: one is the missing mass in 
the magnitude range of the observations, Mi , and another is 
the missing mass in the "whole" magnitude range, Mj. For 
this purpose it is useful to make some considerations on the 
GCS luminosity function. 

F98b reported B and / photometry for ~ 300 globu- 
lar clusters candidates in NGC 1379. This sample was de- 
tected with 4 WFPC2 chips covering a total area of 4.8 
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arcmin . The sample with B < 26 is almost complete; at 
B > 26 the completeness begins to decrease very sharply. . 
The principal source of contamination by foreground stars 
and background galaxies is due to compact, spherical galax- 
ies and F98b provided a background-corrected luminosity 
function as a gaussian function with a peak magnitude of 
(B) = 24.95±0.30 and a width a B = 1.55±0.21. This lumi- 
nosity function is almost complete and uncontaminated in 
the range 21 < B < 25.5. Since we want to estimate the mass 
lost by the GCS, we had to convert the magnitude range 
21 < B < 25.5 in a mass range. This was done assuming 
for the galaxy a distance modulus m — M = 31.32 for NGC 
1379 (Madore, 1986) and , for the GCS, an average colour 
(B — V) = 0.7 and a mass to light ratio (M/L) v& = 1.5 , 
thus obtaining (5.2 x 1O 4 M < m < 3.27 x 1O 6 M ). Adopt- 
ing the mass spectrum (3), the mean mass value of the clus- 
ters (over the told mass rangei) is (m) = 3.5 x 1O 5 M ,. 
This value, together with the number of the missing globu- 
lar clusters (380), leads to Mi = 1.3 x 1O 8 M as estimate of 
the mass fallen to the galactic centre. 
A discussion of the error in Mi induced by the error 
for this agalaxy, as well for NGC 1399 and NGC 1404, 
is postponed to Sect. 3. With regard to the value of the 
mass lost, we remind it represents an understimate of the 
mass removed from the GCS during its evolution because 
it is limited to the actual magnitude range of the obser- 
vations. Hence to compute the "total" mass lost we set a 
minimum and a maximum mass, chosen as m m i n = 1O 4 M 
and m max = 1O 7 M . Actually, normalizing the complete 
mass spectrum to reproduce the number (132) of clusters in 
the observed range (21< B < 25.5), we obtain Ni ttot value of 
179. Moreover, the averaged mass value over the whole mag- 
nitude range is 2.9xlO 5 M thus to give as estimate of the 
total missing mass Mi >to t = 1.5 x lO 8 M ,i.e.l5 % greater 
than the value obtained in the magnitude range of the ob- 
servations. 



2.2 NGC 1399 

NGC 1399 is a giant elliptical galaxy lying about at the 
centre of the Fornax cluster. Among the galaxies around the 
cluster, NGC 1399 is the most populous in globular clusters. 
Its GCS was deeply studied in the past through ground- 
based observations. Among these studies we mention that 
by Kissler-Patig et al. (1997), who deduced a total number 
of globular clusters of 5940 ± 570, value that corresponds to 
the high specific frequency Sn = 11 ± 4. 

In F98a is presented the optical imaging of this galaxy 
obtained with the HST WFPC2 and the Cerro Tololo Inter- 
American Observatory 1.5m telescope. Combining these 
data they provided the spatial distribution of globular clus- 
ters used in the present study. 

2.2.1 The data and their fitting formulas 

The imaging from HST led F98a to set as completeness limit 
B — 26.5. In order to avoid misleading revelations, they 
excluded also all the sources brighter than B = 21, since 
they would be more luminous than known normal globular 
clusters and so, probably, foreground stars. Moreover, these 
authors introduced a colour cut to the data: 
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Figure 2. Surface density profile for the GCS of NGC 1399 (black 
dots) from HST data in F98a. Our modified core model fit is given 
by the solid line, while the dashed line represents the power law 
fit to data given in F98a. The underlying stellar light distribution 
arbitrarly normalized in the vertical direction to match the glob- 
ular cluster system in the outer region is shown as dotted line. 
Radial distance is in arcsec (1 arcsec = 85 pc); surface density in 
arcmin -2 . 



only the sources with 1.2 < B — I < 2. 5 have been in- 
cuded in their sample. We see that 95% of all sources lies 
inside this colour selection. F98a drawn a surface density 
profile for their sample of globular clusters. To do this, they 
introduced a small correction (~ 1.5 %) to take into ac- 
count the number of undetected globular clusters, that are 
expected to be at the faint end of the luminosity function 
(B __; 26.5). Next, they subtracted off the estimated back- 
ground density of 2.9 objects per arcmin 2 . Excluding the 
innermost density points, they fitted the data with the func- 
tion: 

E (r) = 126r _1 ' 2±0 ' 2 (5) 

where r and E(r) are expressed in arcmin and arcmin -2 , re- 
spectively. It is clear from Fig.2 that the GCS surface den- 
sity flattens out in the inner region. The estimate of the 
core radisu given in F98a is 40 arcsec (3.4 kpc). We note 
that this value of the core radius lies in the scatter of the 
relation found by Forbes et al. (1996) between the GCS core 
radius and the galaxy luminosity. 

Assuming 40 arcsec as core radius, our modified-core 
model best fit for NGC 1399 GCS is that with E = 234.42 
arcmin -2 and 7 = 0.61 (solid line in Fig.2). As done first by 
Hanes & Harris (1986) and later by F98a, we assumed 10 
arcmin as limiting radius of the GCS. However it must be 
taken into account that an error on the limiting radius leads 
to an error in the evaluation of the globular cluster number. 
F98a estimated the error in the limiting radius ±1 arcmin; 
consequently, integrating our fit within 10 arcmin, we find 
that the total number of globular clusters in NGC 1399 is 
6100 ± 540. This value agrees well with the values found in 
previous studies (Kissler-Patig et al. 1997 gave 5940 ± 570 
and F98a obtained the value 5700 ±500, adopting the profile 
(6) completed by the inner available data). 
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2.2.2 Number of missing globular clusters in NGC 1399 

The presumed flattening of the GCs distribution in NGC 
1399 has been deeply discussed. While Bridges et al. (1991) 
stated that the GCS has the same slope of the halo stars, 
Wagner et al. (1991) and Kissler-Patig (1997), on the con- 
trary, noticed a flattening of the GCS radial profile toward 
the galactic centre. A weighted average of the slopes pro- 
vided by studies earlier than that of F98a brings to a value 
of -1.47 ± 0.09 for the GCS and -1.72 ± 0.06 for the halo 
stars. F98a found a difference of 0.4 ± 0.22 in the slopes be- 
tween their fitted GCS distribution profile and that of the 
stellar halo, concluding that the flattening of the GCS den- 
sity profile is a real feature. To provide an estimate of 'miss- 
ing' GCs in NGC 1399, we adopted for the stellar light the 
following profile 

log J(r) = -1.6 log r + 2.502, 

obtained shifting the stellar halo density data (Forbes 1999, 
private communication) in the vertical direction, such to 
match the GCS profile in the outer regions. Extrapolating 
inward this profile to the galactic radius of 0.1 arcmin and 
integratig this surface density distribution in the galactic 
region where the GCS and the halo stars profiles differ, i.e. 
~ 0.1 — 8.3 arcmin, we found a number of stars of about 
9680. On the other side, in this radial range we find a num- 
ber of GCs equal to 5165, which yields to an estimate of 4515 
clusters lost, i.e. about 47% of the initial GC population. 



2.2.3 Data completeness and evaluation of the mass lost 

The previous results are limited to the magnitude range 21 < 
B < 26.5 and to the colour range 1.2 < B — I < 2.5. We now 
extend them to fainter magnitudes, obtaining two different 
evaluations of the mass lost in form of GCs. To this scope, it 
is necessary to make some considerations about the colour 
indices of the GCs in NGC 1399. Those in the outer region 
of the galaxy are bluer than those in the central region: 
B - I = 1.70 and 1.97, respectively (F98a). F98a made an 
estimate of the influence of the contamination of background 
galaxies to the globular cluster sample. In this way, they 
found 1.84 as an average value of B — I. This allows us to 
find a value for (B — V), which is needed to transform the 
magnitude interval in a mass range. We have adopted the 
relation (B - V) = 0.347(5 - J) + 0.163 (Grillmair et al. 
1999), obtaining (B-V)~ 0.8. 

Adopting a distance modulus m — M = 31.2 (Jacoby et 
al. 1992) and a colour index (B-V) ~ 0.8 for NGC 1399, and 
a mass- luminosity ratio, for its GCs, equal to 1.5, we obtain 
that the interval 21 < B < 26.5 transforms into the mass 
range 2 x 1O 4 M < m < 3.2 x 10 6 M Q . The mean value of 
the GC mass is then (m) = 2.86 x 1O 5 M . Since the number 
of lost globulars is 4515, the mass loss is 1.3 x 10 9 M Q . 

To estimate the mass lost in the whole GC magnitude 
range, we adopted as total range of globular cluster masses 
that between m m i n = 10 4 M Q and m max = 10 7 Mq. Normal- 
izing the mass distribution (3) to the number of globular 
clusters observed in the range 21 < B < 26.5, we obtain 
a total number of 5686 clusters, to compare with the 5165 
observed in the magnitude range covered by the observa- 
tions. This leads to an estimate of the total mass lost of 



1.44 x 10 9 M Q , which differs of about 10% from that previ- 
ously obtained. 

2.3 NGC 1404 

NGC 1404 is an El galaxy, about a magnitude fainter than 
NGC 1399, sited at about 10 arcmin south east of NGC 
1399. 

As for NGC 1399, much work has been done on the GCS 
of this galaxy. The estimates of the number of globular clus- 
ters of NGC 1404 are made difficult by its vicinity to NGC 
1399. Grillmair et al. (1997) speculated that a great number 
of the NGC 1399 GCs (those of intermediate metallicity) 
have been tidally stripped by NGC 1404. The difficulties in 
distinguishing the two GCSs did not allow to obtain precise 
results from earlier studies, based on photographic plates 
(e.g., Hanes & Harris 1986 and Richtler et al. 1992). This 
problem seems to have been solved by the recent studies by 
F98a, based on data obtained by the HST WFPC2 and the 
1.5 m CTIO telescope. 

2. 3. 1 The data and their fitting formulas 

As for NGC 1399, the magnitude and colour intervals of the 
GC sample are 21 < B < 26.5 and 1.2 < B-I < 2.5. About 
90% of the observed clusters fall in these ranges. Since the 
number of background objects in the field is 14, F98a con- 
cluded that the completeness of the sample of clusters of 
NGC 1404 is > 93%. The GCS surface density profile of 
NGC 1404 has been obtained by subtracting both the back- 
ground density and the density of those clusters which are 
expected to belong to NGC 1399. Excluding the innermost 
data point, F98a fitted the resulting data with the following 
surface density profile: 

E(r) = 36r -1 ' 3±() ' 2 

where r and E (r) are in arcmin and arcmin - , respectively. 
Fig. 3 shows a central flattening of the GCS distribution, 
which allows F98a to define, also in this case, a 'core' for 
the system. Their estimate of the core radius is of about 30 
arcsec (2.5 kpc). As for NGC 1399, this core radius lies in 
the scatter of the correlation found by Forbes et al. (1996). 
We fitted the GCS density distribution data with a modified 
core model, as we have done for the two other galaxies. The 
best fit is that with Eo = 164.39 arcmin -2 , 7 = 0.85 and 
r c = 30 arcsec (solid line in Fig. 3). To obtain an estimate of 
the total number of globular clusters in this galaxy, we have 
integrated our density law in the galactic region within 250 
arcsec (about 4 arcmin), which is the limiting radius of the 
observed data. F98a, in fact, showed that in the direction of 
NGC 1399, the GCS surface density is clearly dominated by 
this galaxy beyond 250 arcsec. 

The modified core model gives a number of globular 
clusters of 744 ±120 in this galactic region (the error in the 
GC number has been calculated adopting an error in the 
limiting radius equal to ±1 arcmin, as suggested by F98a). 
Our estimate of the GC total number is close to that of 
F98a; indeed, using the density profile (8) together with in- 
ner data points they obtained a number of GCs equal to 
725±142. Our value is between that estimated by Harris & 
Hanes (1986), 190±80, and that by Richtler et al. (1992), 
880±120. 
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Figure 3. Surface density profile for GCSs in NGC 1404 (black 
dots) from HST data in F98a. The modified core model fit 
adopted in our study is represented by the solid line, while the 
power law fit to the same data by F98a is shown as the dashed 
line. The underlying stellar light (dotted line) arbitrarly normal- 
ized in the vertical direction to match the globular cluster system 
is also shown. Radial distance is in arcscc (1 arcsec = 85 pc); 
surface density in arcmin - 2 . 



2.3.2 Number of missing clusters in NGC 1404 

Also in this case, to compute the number of globular clus- 
ters lost to the galactic centre of NGC 1404 it is necessary to 
discuss the stellar halo profile. F98a determined a logarith- 
mic slope of the star profile of -1.9 ±0.1. It is clear that the 
NGC 1404 GCS distribution is flatter than that of the halo. 
Thanks to data provided to us in private communication by 
Forbes (1999) we adopted the following halo star profile 



log J(r) = -1.9 log r + 1.76 



(6) 



where r is in arcmin and I(r) in arcmin" 2 , which, as usual, 
was obtained by a vertical normalization so to match the 
GCS profile in the outer regions (see dotted line in Fig. 3). 
Extrapolating this profile inward to the galactic radius of 0.1 
arcmin and integrating this surface density distribution in 
the region of difference with the GCS profile (~ 0.1 — 2.32 
arcmin), we get 1061 as "initial" GC number. Integrating 
our GCS density data fit in the same annulus, we find 508 
presently observed clusters, then we conclude that 553 GCs 
disappeared, i.e. ~ 52% of its initial clusters population. 



2.3.3 Data completeness and evaluation of the number 
and mass of lost globular clusters 

The data analysed are limited to the magnitude range 
21< B < 26.5, transformed into a mass range in order to 
estimate the mass lost from the GCS of NGC 1404. For 
this aim let us do some considerations on the colour in- 
dices of the GCS. According to the analysis of F98a, the 
HST data show that NGC 1404 is dominated by red globu- 
lar clusters, with a significant tail to the blue. Their colour 
distribution yields a mean value {B — I) =1.96, but F98a 
considered necessary to account for the background contam- 
ination. Using the CTIO data they found a bluer colour dis- 



Table 1. The presently observed number of clusters (TV), its ini- 
tial value (iVj), its fractional variation (A), the mass lost in the 
magnitude range of the observations (Mi) and in the whole mag- 
nitude range (Mj ]tot ) 



Galaxy 


N 


Ni 


A 


M,(M Q ) 


M lttot (M Q ) 


NGC 1379 
NGC 1399 
NGC 1404 


132 

5165 

508 


512 

9680 

1061 


0.74 
0.47 
0.52 


1.30 x 10 s 
1.30 x 10 9 
1.37 x 10 s 


1.50 x 10 8 
1.44 x 10 9 
1.75 x 10 8 



tribution, with a mean value (B — I) — 1.71. Now, thanks 
to the relation between B — V and B — I given by Grill- 
mair et al. (1999) we transformed {B — I) into {B — V), 
obtaining a {B — V") value of about 0.75. This value for 
(B — V) together with a distance modulus m — M — 31.2 
maps the magnitude range 21< B < 26.5 onto the mass 
range 1.95 x 10 4 M Q < m < 3 x 1O 6 M . This allows to es- 
timate a GC mean mass (m) = 2.48 x 10 M© and, being 
the number of missing clusters equal to 553, we deduce an 
amount of the mass lost by GCS of 1.37xl0 8 M Q . 

With regard to the total mass lost we have again 
adopted the mass range 1O 4 M < m < W 7 Mq as total 
mass range of globular clusters. Again, normalizing the mass 
spectrum (3) to reproduce the number of globular clusters 
actually observed in the magnitude range 21 < B < 26.5, 
we find 556 as "total" number of GCs. This number, to- 
gether with our estimate of the mean mass value in this 
range (2.9xlO 5 M ), gives a total mass lost of 1.75xl0 8 M Q , 
value 27% greater than the previous one. 



3 THE ERROR ON THE MEAN MASS VALUE 

Using the expressions deduced in Appendix A we can give 
the error in the estimate of the average value of the mass, 
A ^ , of the GCs in the 3 galaxies studied. 



We get 

Aim) Asn Asi As 2 A/3 , s 

_A_/=_ .05 — °_1.48 — --0.21 — 1+0.45-f-, 7 

(m) s si s 2 p 

^) = _ .l5^-1.63^i- 0.19** +0.37^, (8) 

(m) s si s 2 p 

^ = _0.58^- 1.95^-0.19^+0.36^, (9) 

(m) s si s 2 p 



for NGC 1379, NGC 1399 and NGC 1404, respectively. From 
these expressions it is clear that the greatest source of inde- 
termination is due to si, i.e. to the parameter of the mass 
spectrum that characterizes the intermediate mass range. 
However it is important to underline that the abundance of 
the GCs sample in these ranges allows an accurate determi- 
nation of si. Assuming an error of 10% in the parameters 
so, S2 and f3 = (M/L) v& , we obtain, as resulting errors in 
(m) 22%, 23% and 30%, for NGC 1379, NGC 1399 and NGC 
1404, respectively. 
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4 THE CORRELATION BETWEEN THE MASS 
LOST BY GCS AND THE CENTRAL 
GALACTIC BLACK HOLE MASS 

Even if theoretical arguments and observational data arc 
clearly compatible with the explanation of the difference 
between the GCS and stellar bulge radial distributions as 
a result of dynamical evolution, we still cannot firmly state 
it. By the way it is interesting, and surely constitutes an 
argument in favour of this interpretation, to see how our 
results for the (supposed) mass loss in form of glbobular 
cluster correlates with the galactic central black hole mass, 
at least for the set of seven galaxies for which these data are 
available (our galaxy, M 31, M 87, NGC 1399, NGC 4365, 
NGC 4589, IC 1459). 

Sources of data in Fig. 4 are: for Mi this paper and Capuzzo- 
Dolcetta & Vignola, 1997, while the black hole masses are 
from Magorrian et al. (1998) except for our Galaxy and IC 
1459 taken from Gebhardt et al. (2000). A recent estimate 
by Saglia et al. (2000) gives as compatible with kinematic 
data of the inner 5 arcsec of NGC 1399 a black hole mass 
~ 5 x 10 8 Mq, i.e. about a factor of 10 less than the value 
estimated by Magorrian et al. (1998) and reported in Fig. 4. 

Fig. 4 shows, even in the paucity of data available, a 
clear increasing trend of Mi as function of Mbh- Actually, 
the least square straight line fitting the data relative to 
galaxies with black hole masses heavier than 10 s M© has 
slope 0.99, i.e. the Mi - Mbh relation is remarkably linear, 
while the straight line connecting the 2 low mass points has 
slope 0.08. The best simple fit to the whole set of data is 
given by LogMi = a + bexp(LogMth) with a = 7.158 and 
b = 1.194x 10~ 4 , giving r 2 = 0.76 and standard error = 0.44 
(the least square straight line fitting the whole data set 
yields r 2 = 0.61 and standard error = 0.56). Incidentally, 
we note that the steeper slope of the Mi - Mbh relation for 
high values of the nucleus mass is what expected when the 
black hole is very massive and it acts as powerful engine of 
GCS evolution. 



5 CONCLUSIONS 

We used data taken with the WFPC2 of the Hubble Space 
Telescope to compare the globular cluster system distribu- 
tion to that of bulge-halo stars in a sample of three elliptical 
galaxies in the Fornax cluster. 

The first relevant outcome of this study is that the ra- 
dial distributions of the globular cluster systems in NGC 
1379, NGC 1399 and NGC 1404 confirm that the GCSs 
in elliptical galaxies are less concentrated to the inner re- 
gions than the halo stars. Following the scheme of previous 
works (McLaughlin 1995, Capuzzo-Dolcetta & Vignola 1997, 
Capuzzo-Dolcetta & Tesseri 1999), i.e. assuming that the 
globular clusters and halo star distributions were initially 
the same and the present difference is due to evolutionary 
effects acting on the GCS, we computed the number of glob- 
ular clusters lost as the integral of the difference between the 
two surface distributions (suitably normalized). 
In this way we have found that the GCSs of NGC 1379, 
NGC 1399 and NGC 1404 could have lost the 74%, 47% 
and 53% of their initial populations, respectively. These re- 
sults are, obviously, limited to the magnitude range of the 
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Figure 4. The correlation between the (logarithmic) GCS mass 
lost and the central galactic black hole mass for 7 galaxies for 
which these data are available. Masses are in solar masses. The 
straight solid line is the least square fit to all data, the dashed 
line is the least square fit restricted to data of the most massive 
black holes (LogMbh > 8) while the curve is the exponential fit 
described in the text 



observations within which the data reach a sufficiently high 
level of completeness. 

The large number of 'missing' globular clusters obtained re- 
veals that a significant amount of mass could have reached 
the galactic centre, too. An estimate of the mass fallen to 
the galactic center needs the assumption of a value for the 
initial mean mass of globular clusters, (m) . It is reliable (see 
Capuzzo-Dolcetta & Tesseri 1997) to assume for (m) the 
mean mass value of the presently observed globular clus- 
ters. This quantity can be estimated assuming a mass spec- 
trum, that we took as a piece-wise power law for it seems 
to represent one of the best approximations for elliptical 
galaxies globular cluster systems (McLaughlin 1994; Har- 
ris & Pudritz 1994). In particular, we adopted as unique 
mass spectrum for the three galaxies of the Fornax cluster 
that obtained by normalization of the luminosity function 
adopted by McLaughlin (1995) for the Virgo elliptical giant 
M87, and assuming a mass to light ratio (M/L) v& =1.5. 

We provided two different estimates of the mass lost in 
form of globular clusters: the first is that due to clusters in 
the magnitude range of the observations and the other refers 
the "presumed" whole magnitude range. The first estimate, 
Mi , is simply obtained as the product of the mean mass and 
the number, AT;, of the missing globular clusters. With re- 
gard to the second estimate we made the, quite reasonable, 
assumption that the total number of missing globular clus- 
ters scales with the total number of the clusters present in 
the galaxy. Hence we computed this total number of globu- 
lar clusters assuming 10 4 M Q < m < 1O 7 M as whole mass 
range, so to include the contribution to iV; coming from un- 
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observed clusters. The lower limit of the mass range has not 
a big influence on the total mass; with regard to the upper 
limit, the large decreasing slope (—3) of the MF at large 
masses ensures a high level of convergence of the total mass 
reached when assuming 10 7 Mq as maximum mass. The val- 
ues of the total mass lost, Mi, obtained in this way for the 
three galaxies NGC 1379, NGC 1399 and NGC 1404 are 
not too different from those estimated in the mass ranges 
of the observations (the difference is about 15 %, 10 % and 
27 %, respectively; see Table 1). This is explained by that 
the range of magnitude covered by the observations is close 
enough to the whole magnitude range, thanks to the high 
power of HST. 

It is relevant noting the positive correlation we found, 
for the available set of 7 galaxies in Fig. 4, between our 
evaluated mass lost from the GCS system and the central 
black hole mass of the parent galaxy: this is what expected 
on the basis of globular cluster-black hole tidal interaction 
because the more massive the central black hole the greater 
its influence on the GCS population. 

To conclude, we remark how the values of the mass 
lost in form of decayed and destroyed globular clusters is 
significant and so it should have been played a role in the 
central region activity of their parent galaxies. 
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APPENDIX A: 

A significant source of error on the evaluation of the mass 
removed from the GCS is due to the assumption of the ex- 
ponent s of the mass spectrum and of the value of the mass- 
light ratio (M/L) e (hereafter 0). Infact, an error on these 
parameters leads to the relative error on (m) 



A(m) _ d{m) s As d{m) (3 A/3 
(m> 



ds (m) s ^ dj3 (m) (3 



(Al) 



Before giving an explicit expression of (Al), we note that, if 
we use the mass spectrum as in (3) , (m) is so defined: 



r max k(L)L' s+1 dL 

(m) = (PL) = P(L) = p JL — 

J L max k(L)L-"dL 



(A2) 
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that is, in terms of mass 

Im k(m)m~ s+1 dm 



(m) 



XT" 10 " fc( m ) m ~ s ^ m 



(A3) 



where k(m) is obtained by the proper tranformation of the 
McLaughlin (1994) normalizing factor k(L), i.e. 



k(m) 



A 

Ami 1 - 30 
Am ei - S0 m s 2 2 - S1 



m m in < m < mi 

mi < m < mi . (A4) 

m 2 < m < m max 



where, of course, m m i„, m max , mi and m,2 depend linearly 
on p. To give (Al) a more compact form, we have put 



m = m min and m 3 = m max , 
fco = A, ki = Ami 1 ' 30 , k 2 = Ami 



Ntot = f "' max k(m)m s dm. 

Jm min V ' 

Remembering that s — (so,si,S2) (Al) becomes 



d(m) s, As, _^ d(m) p Ap 
'"" ^ dp (m) P 



(m) ^ dsi (m) Si 

i= 

where 9 ^ , are so defined: 



d(m) 
dsi 



—y 

Ntot ^ 



3 = 



(j + 1,2)-^ (j,2)) + 
-<m) (tf 0" + 1,1)- ^ (j,l)) 



(A5) 



(A6) 



with 

if \ 9kj m„ 



+ S 3 - — — : — kj ( — —. In m r ) ,and 



2 r ( 

d(m) 1 \ ^ kj y 



9P 



m j+ + -m 3 



'Vl- 1 -'V v '" 



.(A7) 
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